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Abstract 

Photoproduction data collected with the HI detector at HERA in 1994 are used to 
study the cross-sections for inclusive charged particle production and the structure 
of the photon. The differential cross-sections da/dp T , for \rj\ < 1 in the HERA labo- 
ratory frame, and da/dr] for px > 2 GeV/c and px > 3 GeV/c have been measured. 
Model calculations of these cross-sections, based on perturbative QCD, indicate that 
the results are sensitive to the parton densities of the photon as well as to higher 
order effects, which are phenomenologically treated by multiple interactions. This 
sensitivity is exploited to determine the leading order x 7 distribution of partons in 
the photon using a new method based on high px charged particles. The gluon 
content of the photon is extracted and found to rise with decreasing x 7 . 
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1 Introduction 



The observation that the photon has a hadronic structure was first made more than 30 years 
ago when the cross-sections of hadronic photoproduction interactions were demonstrated to have 
dependences on energy and momentum transfer which were similar to those in hadron-hadron 
interactions Q. With the advent of the quark-parton model, and subsequently QCD, more 
quantitative predictions for this hadronic structure became available, and its gross features were 
identified experimentally in e + e~ interactions Subsequently these features were also observed 
in hard photoproduction processes [3-15]. 

Measurements of the photon structure functions in e + e~ interactions are directly sensitive 
to the quark structure of the photon. Only through QCD evolution studies can information 
be extracted concerning the gluon component of this structure, but the presently available 
data have not been precise enough for such an analysis. Recently, studies of jets and high px 
charged particles in photoproduction events at the ep collider HERA have shown sensitivity to 
the partonic content of the photon. Here the photon structure is probed by the partons of the 
proton, rather than by a virtual photon as in ej collisions. Hence these data are sensitive to both 
the quark and gluon content of the photon. Leading Order (LO) diagrams are shown in Fig. [l], for 
so called direct (Fig. la; the photon couples directly to the partons in the proton) and resolved 
(Fig. lb; the partons of the hadronic component of the photon scatter on the partons of the 
proton) processes in photoproduction. In LO QCD only the latter process contains information 
on the partonic structure of the photon. 

In H inclusive charged particle production has been studied. It was established that the 
data, in particular the tail at large transverse momentum, pr, can be described by Next-to- 
Leading Order (NLO) QCD calculations. The charged particle distributions as a function of the 
pseudorapidityQ, rj = —In tan(0/2), were found to be sensitive to the partonic structure of the 
photon. In this paper we present differential cross-sections dajdp\ for \rj\ < 1 in the HERA 
laboratory system, and da/drj for px > 2 GeV/c and pt > 3 GeV/c for charged particles, 
measured with the HI detector. Photoproduction events have been selected by tagging the 
scattered electron. The data are based on an event sample which is 50 times larger than the one 
used in j5j. 

In a first measurement of the gluon content of the photon was presented, while in || 
an effective parton density for the photon was extracted. In this paper the distribution of the 
momentum fraction x 7 = E parton / E 1 of the parton of the photon entering the hard scattering 
process shown in Fig. lb is measured and used to extract the gluon density in leading order 
in the photon. This analysis follows closely that presented in Q, where the kinematics of the 
hard scattering process were reconstructed by means of jets measured in the detector. In the 
analysis presented here the measured variable x^ ec , which is found to be correlated with the 
true value of x 7 of the parton entering the scattering, is based on charged tracks with a high 
transverse momentum pt- This method thus avoids two large systematic errors entering in the 
jet analysis. These are the energy scale uncertainty of the calorimeter and the uncertainty of 
the jet energy measurement due to overlap with energy deposits from soft multiple interactions 
which may occur on top of the hard scattering process. The drawback of this method is a 
stronger sensitivity to fragmentation uncertainties. 



1 9 is the polar angle of the particle in the HERA laboratory frame, measured with respect to the proton beam 
direction. 
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Figure 1: Examples of diagrams for direct (a) and resolved (b) photon processes in ep scattering; 
in (b) a gluon has the space-like virtuality of the probe and a quark is resolved in the photon; 
similarly gluons in the photon can be resolved by quarks and gluons with space-like virtuality. 

2 The Apparatus 

The data were collected with the multi-purpose detector HI at the HERA collider in 1994, in 
which electrons^] of 27.5 GeV collided with protons of 820 GeV. The total data sample corre- 
sponds to a luminosity of 1.35 ± 0.02 pb _1 . 



A detailed description of the Hl-detector has been given elsewhere [16]. Here only the 
components crucial for this particular analysis will be briefly described. 

The HI luminosity system consists of an electron tagger and a photon tagger, located 33 m 
and 103 m from the interaction point in the electron beam direction, respectively. The lumi- 
nosity is determined from the rate of Bremsstrahlung process ep — * epj events by detecting the 
photon. Both detectors are TICl/TlBr crystal Cerenkov calorimeters with an energy resolution 
of 10%/y E(GeV). The electron detector allows to tag the photoproduction events used in this 
analysis by detecting electrons scattered at small angles. 

The central tracker (CT) consists of inner and outer cylindrical jet chambers, z-drift cham- 
bers and proportional chambers. The two cylindrical drift chambers | ]17| (CJC), are mounted 
concentrically around the beamline inside a homogeneous magnetic field of 1.15 Tesla and pro- 
vide up to 65 space points for a charged track, yielding particle charge and momentum from 
the track curvature in the polar angular range of 20° and 160°. The transverse momentum 
resolution is ct(pt)/pt ~ 0.6% -pr (GeV/c) in the px range 2 < pt < 12 GeV/c considered here. 
The resolution for the polar angle 9 is 2 mrad. 

The central tracking system is complemented by a forward tracking system. All trackers are 
surrounded by a fine grained liquid argon sampling calorimeter, consisting of an electromagnetic 
section with lead absorbers and a hadronic section with steel absorbers. 



3 Data Analysis 



2 In 1994 the incident lepton at HERA was a positron, but we keep the generic name "electron" for the incident 
and scattered lepton throughout this paper 
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3.1 Event Selection 

Photoproduction events were selected which have the scattered electron detected in the elec- 
tron tagger, and consequently have a four-momentum transfer Q 2 < 10~ 2 GeV 2 , where Q 2 = 
AE e E' e cos 2 (9/2). Here E e and E' e are the incoming and scattered electron energy and 9 the polar 
angle of the scattered electron. Tagged events were required to have an electron candidate in 
the fiducial volume of the electron-tagger with energy E' e > 4 GeV and to have less than 2 GeV 
deposited in the photon detector. The latter condition suppresses background from the proton 
beam, which appears in a random coincidence with the high rate Bethe-Heitler events, and also 
reduces QED corrections. In addition a track-based trigger was required, which is formed using 
only 10 selected layers out of 56 radial signal wire layers of the CJC. It demands the presence 
of at least one track of negative charge with px > 700 MeV/c and with a distance of closest 
approach of less than 2 cm from the nominal beam axis; this requirement suppresses beam-wall 
background as well as fake tracks from random hits due to synchrotron radiation. 

For this analysis the events were required to have 0.3 < y = E^/E e < 0.7 where E 7 is the 
photon energy. As a consequence, in the jp centre-of-mass system (CMS) the energy range 
is 165 < ^/s 7P < 251 GeV, with an average of ^JSryp ~ 200 GeV. An interaction vertex, 
reconstructed using tracks in the CT, had to be in the range —25 < zy < +35 cm, where zy is 
the z coordinate of the reconstructed vertex. 

A total of about 960 000 events satisfy all described criteria. The remaining non-ep back- 
ground was determined by analysing data from electron pilot bunches (i.e. electron bunches 
which have no colliding proton bunch partner) and from a monitoring trigger of minimum bias 
events. It was found to be less than 2.7% at this stage of selection, before tighter cuts on track 
quality and transverse momentum were applied. 

3.2 Track Selection 

For this analysis tracks from the central tracker were selected which have a transverse momentum 
with respect to the beam axis of at least 2 GeV/c, a pseudorapidity measured in the HERA 
laboratory reference system of \r]\ < 1 and a minimum track length AR of 30 cm in the x — y 
plane. The latter cut is essential to ensure a good measurement of pt of the track, from 2 up 
to 12 GeV/c, which is the range considered here. To minimize the systematic errors, sectors 
of the CT in (f>, the azimuthal angle, which had a lower efficiency have been excluded from the 
analysis. A total of 16591 tracks was selected for the analysis, resulting from 15543 events. 

3.3 Corrections and Background 

The residual background from beam-gas, beam-wall, and pile-up events, after demanding a well- 
measured, high pt track such as defined above, was found to be negligible. The background 
from cosmic events was removed by rejecting events with two oppositely charged tracks meeting 
the requirements above (except for region cut), and which were aligned in 9 and <f> to better 
than 5°. All events rejected by this cut were visually checked and confirmed to be cosmic events, 
and all were found to have at least one track with pt > 6 GeV/c. Conversely, all events with 
at least one track with px > 6 GeV/c were visually scanned but no further cosmic event was 
found. A check for the potential contamination from DIS overlap events was also performed: no 
such events were found. 

In order to obtain the produced number of tracks, N tr , the observed number of events and 
tracks were corrected as follows: 
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• A etag : The geometrical acceptance of the electron tagger, determined as described in IS], 
was corrected as a function of y. On average the acceptance is 55% for 0.3 < y < 0.7. The 
tagger efficiency within the chosen acceptance is 100%. 



• etrig- The CJC trigger efficiency has been estimated by comparing the detection rates for 
samples of events with the main trigger and a minimum bias trigger for e-tagged events. 
For the high pr range considered here, it was estimated to be 96.3 ± 1.2%, independent of 
Pt, V and cf). 

• etrack '■ The overall track efficiency contains three multiplicative factors : (i) The efficiency 
for finding and reconstructing single tracks, e rec , was taken from Monte Carlo studies and 
found to be 99.2 ± 0.8%. The agreement between data and Monte Carlo was verified by 
visual scanning of events in both samples, (ii) The efficiency of the cut on transverse 
track length, £a_r, was estimated from data to be 97.9 ± 0.4%. This was checked to be 
independent of both 7] and (p. (iii) The (^-restriction, e^, of the geometrical acceptance 
of the central drift chamber was accounted for assuming a uniform distribution in (j>, and 
amounts to 53.1%. 

• The bin-to-bin migration in pTi tmigi due to the rapidly decreasing pt spectrum and to 
the finite resolution in p?, was estimated from a Monte Carlo simulation and for the 
varying bin sizes used here it was found to be a simple multiplicative factor 103.8 ± 1.3%, 
independent of pr- 

The total correction factor for tracks is then 51.5 ± 2.0%, where the systematic errors 
have been added in quadrature. The uncertainties on these correction factors, were taken as 
systematic errors. Detailed checks and comparisons with Monte Carlo led to the introduction 
of an extra systematic error, to account for a remaining uncertainty in the pt measurement, 
which showed up in two ways: in the ratio of positively and negatively charged tracks and in 
the efficiency determination for a stronger cut on the transverse track length (AR > 50 cm). 
The resulting uncertainty was parameterized as 2% • ^/(GeV/c) at the centre of each bin for 
Pt distributions, and as 2% ■ p 1 ™ /(GeV /c) for i] distributions, with p 1 ™ the lowest pt value of 
particles included in the distribution. 

All systematic uncertainties were added in quadrature, except for the overall uncertainty of 
5% from the luminosity measurement and the electron tagger acceptance which is not included 
in the error bars shown in the figures below. 



4 Monte Carlo Simulation Programs 



The PYTHIA 5.7 event generator p| was used to simulate photon-proton interactions. As 
well as the leading order cross-section, PYTHIA includes initial- and final-state QCD parton 
radiation effects which are calculated in the leading logarithmic approximation. The strong 
coupling constant a s was calculated in first order QCD using Aqcd = 200 MeV for 4 quark 



flavours. For the parton distributions of the proton the GRV-LO [2C] set was used. 



The model calculations were made for different sets of the parton densities in the photon: 



GRV [J21], SAS1D |22| and LAC1 |2J. All these parametrizations were extracted from QCD 
fits to photon structure function data measured in collisions from e + e~ interactions, but 
have different assumptions for the gluon content, which is only weakly constrained by these 
measurements. In particular LAC1 assumes a large gluon content of the photon at small- x 
values which is larger than for both other parametrizations. The GRV and SAS1D distributions 
both start the evolution from a small starting scale, Qq = 0.25 GeV 2 and 0.36 GeV 2 respectively, 
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a procedure which has turned out to be quite successful for the parton densities in the proton. 
The different treatment of the vector meson valence quark distributions leads to a larger gluon 
component at small- x of the photon for the GRV compared to the SAS1D parton densities. 

Multiple parton interactions were generated in addition to the primary parton-parton scat- 
tering. They are calculated as leading order QCD processes between partons from the photon 
and proton remnants. The transverse momentum of all (primary and multiple) parton-parton 
interactions was required to be above a cut-off value of j>™ m , depending on the photon structure 
function. For the parton densities used here the values p™ n = 1.2, 1.0 and 2.0 GeV/c have been 
used for GRV, SAS1D and LAC1 respectively. These values have been found to give an optimal 
description of the transverse energy flow outside jets @. 



The PHOJET 1.06 event generator [24|, based on the two-component Dual Parton model 25] 
has been used as well. PHOJET incorporates very detailed simulations of both multiple soft 
and hard parton interactions on the basis of a unitarization scheme. It also includes initial- and 
final-state interactions. For the distributions shown below, it was verified that the PYTHIA and 
PHOJET predictions for the GRV parton distributions for the proton and the photon agree to 
better than 5%. 

Hadronization in both PHOJET and PYTHIA was modelled with the LUND string frag- 
mentation scheme (JETSET 7.4 [p6[ ). Since the kinematics of the hard scattering process will 
be related to the pt values of the particles with high transverse momentum, it is imperative to 
study the fragmentation dependence of this measurement. A model with a different hadroniza- 
tion scheme, HERWIG [^], was used for this purpose. In this model the hadronization is based 
on cluster fragmentation. Other parameters, such as the parton densities for the photon and 
proton, were taken to be the same as for PYTHIA. Both the LUND model and HERWIG were 
found to give a good description of the general features of fragmentation as measured in e + e~ 
collisions at LEP (e.g. 



5 Inclusive Charged Particle Cross Sections 

The invariant cross-section for single particle production is given by 

d 2 a f d 3 cr , , „ d 3 <7 



dp T di] J dp T dr] dcj) r dp 3 ^ ^ 

assuming azimuthal symmetry of the cross-section allowing integration over <p. The measurement 
was made at a .yj 's^p rN ~' 200 GeV, and. effectively averages over the region 165 GeV <C ^ ' s^p <C 
251 GeV. No significant ^/s 7P dependence was found when the data were subdivided in two ^/ s<yp 
bins. The cross-section for inclusive charged particle production in 'yp collisions was calculated 
from the corrected number of tracks produced, Nt r , in a bin of pt and r/. It is given by: 
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Ntr(PT,v) 



(2) 

dp T di] C-F- 2 PT Ap T -Ar) w 

where C denotes the integrated luminosity, F is the photon flux integral and the flux f(y) is 
defined according to dcr(ep) = a(^p) ■ f(y) ■ dy. For the chosen y-range the integral over y of 



the photon flux yields F = 0.0136, assuming the Weizsacker- Williams approximation [29, pCfl . 
At/ and Apj, = 2 ■ px ■ Apx are the bin widths. 

The resulting differential cross-section for the sum of positive and negative charged particles is 
shown in Fig. [2|a. The error bars denote the statistical and systematic errors added in quadrature. 
The measured cross-sections are listed in Table ffl. The data exhibit a strong improvement in 
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Figure 2: a) The inclusive 7p cross-section for charged particles in photoproduction (full 
squares) measured in the kinematical region | rj |< 1.0 at an average -^/s-yp ~ 200 GeV\ The error 
bars denote the statistical and systematic errors added in quadrature. An overall uncertainty 
of 5% from the luminosity measurement and the electron tagger acceptance is not included in 
the errors. The curves indicate power-law fits, as described in the text, for these data and 
for the pp data from UAl and CDF as given in |3^, [33|. b) The ratio of data over the NLO 
QCD calculation with scale for £ = 1. The shaded band shows the expected variation of 
the ratio as £ changes from 0.5 to 2, illustrating the sensitivity of the QCD calculation to the 
renormalization and factorization scales (see text). 
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Figure 3: The differential ■jp cross-section da/dr/ for inclusive production of high px charged 
particles (full squares) in comparison with LO QCD calculation by PYTHIA (histograms). The 
error bars denote the statistical and systematic errors added in quadrature. An overall uncer- 
tainty of 5% from the luminosity measurement and electron tagger acceptance is not included. 
Different lines represent different photon structure function parametrizations: GRV with (full) 
and without (dash-dotted) multiple interactions, SAS1D (dashed) and LAC1 (dotted). The 
contributions from the direct photon processes are shown as shaded histograms. 

precision compared to the previous HI measurement || and cross-sections are now given up to 
12 GeV/c in px- The QCD inspired power-law [31] expression 



E.f^ = A.(l + » (3) 
dp 6 {pt)o 

was fitted to the data. The fitted curve is shown in Fig. ||a, and describes the data well over 
the whole px range. The fit gives A = 5.44 ± 0.66 mb and n = 7.03 ± 0.07 (statistical errors). 
The parameter (px)o was fixed to the value 0.63, as found in ||, but other values in the range 
0.5-0.75 were found to give equally good fits. The effect of the choice of (px)o leads to an 
uncertainty in the power n of about 0.2. 

Also shown in Fig. 0a are the results of similar fits to the pp collider data of the UA1- 



collaboration [32] in the rapidity region | rj |< 2.5, at CMS energies yfs = 200 and 900 GeV, and 
of CDF [33 1 in the rapidity region | rj |< 1, at ^/s = 1800 GeV, scaled by the factors indicated 
in the figure, which essentially normalizes all data to the 7p cross-section at pt = 2 GeV/c. 
The high px tail is observed to increase with increasing energy in pp. The hardness of the "fp 
spectrum is comparable to that of the pp data at the highest energy. The high px tail in the data 
is clearly larger than in pp collisions at similar CMS energy, which can be understood as being 
due to extra contributions in jp, namely the direct and the pointlike resolved (or anomalous) 
component ||. 

In Fig. [2]b the ratio of data to the NLO calculation [34] using GRV structure functions for 
photon and proton is shown. Both renormalization and factorization scales were chosen to be 
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Vt fGeV/cl 

11 L / J 


4^- (nb/(GeV/c) 2 ) 


A-tt-t- (nb/(GeV/c) 2 ) 


2.00 - 2.08 


685.0 


35.2 


2.08 - 2.16 


539.1 


29.3 


2.16 - 2.24 


434.2 


24.7 


2.24 - 2.34 


359.9 


20.7 


2.34 - 2.44 


274.3 


16.7 


2.44 - 2.54 


216.6 


13.9 


2.54 - 2.66 


184.8 


12.1 


2.66 - 2.78 


135.9 


9.53 


2.78 - 2.90 


102.8 


7.64 


2.90 - 3.06 


78.32 


5.97 


3.06 - 3.22 


56.41 


4.59 


3.22 - 3.38 


45.62 


3.97 


3.38 - 3.58 


31.01 


2.82 


3.58 - 3.78 


23.47 


2.32 


3.78 - 4.02 


14.33 


1.49 


4.02 - 4.30 


11.79 


1.33 


4.30 - 4.58 


6.905 


0.885 


4.58 - 4.92 


5.210 


0.692 


4.92 - 5.32 


3.224 


0.465 


5.32 - 5.72 


1.673 


0.283 


5.72 - 6.20 


1.148 


0.225 


6.20 - 6.90 


0.631 


0.126 


6.90 - 7.80 


0.408 


0.100 


7.80 - 8.80 


0.232 


0.068 


8.80 - 10.0 


0.066 


0.020 


10.0 - 12.0 


0.019 


0.007 



Table 1: Measured differential cross-section and total errors (A) for the production of charged 
particles in the r\ range —1 to 1. An overall uncertainty of 5% from the luminosity and electron 
tagger acceptance is not included. 





p T > 2.0 GeV/c 


p T > 3.0 GeV/c 


V 


%w 


^ (*>) 


% ( nb ) 


±% W 


-1.0 - -0.8 
-0.8 - -0.6 
-0.6 - -0.4 
-0.4 - -0.2 
-0.2 - -0.0 
0.0 - 0.2 
0.2 - 0.4 
0.4 - 0.6 
0.6 - 0.8 
0.8 - 1.0 


1431.1 
1473.1 
1556.3 
1593.3 
1611.9 
1643.8 
1730.2 
1721.3 
1790.4 
1721.2 


77.8 
79.4 
83.2 
85.1 
85.3 
87.5 
91.2 
90.9 
94.0 
90.6 


293.9 
312.3 
327.6 
350.3 
337.4 
343.7 
396.8 
346.1 
353.3 
351.4 


27.2 
28.4 
29.3 
30.7 
29.6 
29.8 
34.8 
31.4 
31.3 
31.5 



Table 2: Measured differential cross-section and total errors (A) for the production of charged 
particles. An overall uncertainty of 5% from the luminosity and electron tagger acceptance is 
not included. 
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Figure 4: The differential jp cross-section da/drj for inclusive production of the high pt charged 
particles (full squares) in comparison with NLO QCD calculations |34[] , with different values of 
the scale, taken as fi = ip\. The error bars denote the statistical and systematic errors added in 
quadrature. An overall uncertainty of 5% from the luminosity measurement and electron tagger 
acceptance is not included. 



p\. The effect of this choice is shown by predictions for a scale of with £ = 0.5 and 2, shown 
as a ratio to the prediction for £ = 1. The calculations describe the data well, in particular when 
£ is close to one. 

The dependence of the cross-section on the pseudorapidity rj, da/dr/, is shown in Fig. || for all 
particles with px larger than 2 and 3 GeV/c. The measured cross-sections are listed in Table ^. 
Note, that the rj measured in the HERA laboratory frame is on average shifted by +2 units 
with respect to the rj* distribution in the 7p CMS system. The spread of this shift due to the 
photon energy range is less than 0.3 units in pseudorapidity. Therefore the rj and rf distributions 
look very similar. The data are consistent with being flat in most of the region, with a slight 
decrease towards rj = — 1, i.e. in the photon direction. The data are compared with predictions 
of LO QCD calculations made with PYTHIA, for different structure functions of the photon 
and assumptions on the soft interactions in the underlying event (cf. section 4) . Predictions are 
shown for GRV with (full) and without (dash-dotted) multiple interactions, SAS1D (dashed) and 
LAC1 (dotted). For the latter two predictions multiple interactions were included. All structure 
functions except LAC1 show a falling distribution with increasing ri, contrary to the data. The 
contributions from the direct photon processes are shown as shaded histogram and decrease 
towards the proton direction. For the data with pr > 2 GeV/c, there is a strong sensitivity 
to the parton distributions in the photon, and also to the amount of underlying interactions as 
shown by the GRV curves. None of the predictions presented shows satisfactory agreement with 
the data. For data with pr > 3 GeV/c, the effect of multiple interactions becomes very small, 
as shown for the GRV predictions. The sensitivity to the photon structure is also reduced but is 
still significant. For large pt the data agree best with LAC1, and disagree with SAS1D. It was 
checked that these results do not depend on the choice of the parton distributions of the proton, 
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when selected from those consistent with the most recent proton structure function data. 

In Fig. |] the cross-sections as a function of 77 are compared with NLO calculations, including 
direct and resolved contributions. The NLO calculations J34| reproduce the data, however one 
can see a significant effect of the choice of the QCD scale on both the shape and normalization 
of the distributions, as shown by the curves for different £ values. The effect is smaller for 
data with px > 3 GeV/c. The scale has also been found to be the reason for the considerable 
difference between the LO and NLO GRV cross-section predictions. PYTHIA uses as a scale 
the pt of the hard partonic interaction, while the program of |34}| uses the px of the final state 
particle, which can differ substantially from that of the original parton. 



6 The x 1 Distribution and Gluon Density in the Photon 

In this section the distribution of the momentum fraction x 7 = E var t on j of the parton of the 
photon entering the hard scattering process will be determined from the charged tracks. For this 
analysis events were kept which have at least one track reconstructed with a pr > 2.6 GeV/c. 
This value is a compromise between large systematic uncertainties from multiple interactions 
and the statistical precision. In total 9378 events have been selected. 

For each event with at least one measured track of a charged particle with px > 2.6 GeV/c 
the variable x^ ec was calculated: 

x r ec = ZEE£l (4 ) 

where the sum runs over all tracks with px > 2 GeV/c. This variable is correlated to the true x 7 
of the interaction. Apart from the strongly reduced sensitivity to effects of multiple interactions, 
the requirement px > 2.6 GeV/c ensures that the data stay safely away from the region which 
is affected by the pV^ in cut used in the Monte Carlo event generation (cf. section 4). 

Monte Carlo events were used to study the correlation between the measured x^ ec and the true 
x~, shown in Fig. |5[ The generated events were fed into the detailed HI simulation program and 
then subjected to the same reconstruction and analysis chain as the real data. The correlation 
is shown for the Monte Carlo events generated with the PYTHIA program, using the LAC1 
parametrization of the structure function of the photon. 

For the reconstruction of the true £ 7 an unfolding procedure was used |p5[l , following the 
analysis as detailed in Q. Monte Carlo events generated with PYTHIA using the LAC1 parton 
densities have been used for the unfolding. The result is presented in three bins 111 %ry 3 EL condition 
imposed by the unfolding procedure, in order to minimize the bin-to-bin correlations. The 
resulting x^ distribution is shown in Fig. |(| with statistical errors only, which include the bin-to- 
bin correlations from the unfolding. The dotted curve is the prediction of the PYTHIA model 
with the LAC1 parton distributions for the photon. The dashed curve shows the component of 
the event rate where only quarks from the photon side are involved. The quark distribution in 
the photon is rather well constrained from e^f measurements as discussed in section 1 above. The 
full calculation gives a fair description of the measurement, and a large contribution of gluon 
induced processes (from the photon side) || is clearly confirmed. 

The gluon distribution in the photon, extracted at leading order, was then obtained from 
the x 7 distribution as follows. The contribution of processes in which quarks are resolved in 
the photon (for example Fig. |l|(b)) were calculated and subtracted from the data. The ratio 
of the subtracted data distribution with the distribution calculated from Monte Carlo using 
only processes initiated by a gluon from the photon side gives weight factors. Applying these 
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Figure 5: Correlation between reconstructed and true x 7 values using the PYTHIA Monte 
Carlo with the LAC1 parton densities for the photon. 



weight factors to the input gluon distribution used in the Monte Carlo yields the measured gluon 
distribution, which is shown in Fig. f?|. The average transverse momentum squared of the hard 
partonic scattering for this data sample amounts to < ptf, >= 38 GeV 2 /c 2 , as derived from Monte 
Carlo studies, and is taken as the scale for comparisons with parton densities. The total errors 
include statistical and systematic errors where all contributions have been added in quadrature. 
Apart from the systematic errors already included in the cross-section measurements discussed 
previously, the following systematic errors specific to the unfolding had to be taken into account: 

• The uncertainty in the quark density of the photon was found to be 15% by repeating the 
analysis with different parametrizations of the quark distributions [pl| , 22, 23 1. 



A variation of the unfolding parameters led to a 10% uncertainty. 

Changes of up to 20% were observed when different fragmentation models (PYTHIA, 
PHOJET, HERWIG) and different parton densities for the photon structure (GRV, LAC1) 
were selected in the Monte Carlo used for the unfolding. Hence a 20% systematic error 
from this source was assumed. 



The uncertainty in the parton densities of the proton is negligible compared to the uncertainties 
above and this was neglected. 
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Figure 6: The x 7 distribution in the photon. The errors are statistical only. The curves are 
full (dotted line) and 'quark only' (dashed line) predictions of PYTHIA, using the LAC1 parton 
distributions of the photon. 

The result is compared with the measurements based on the 1993 HI data |J, using jets 
instead of charged tracks. The average transverse momentum squared of the hard partonic 
scattering for the jet data sample amounts to < p\ >= 75 GeV 2 /c 2 . The measurements are 
found to be consistent. The improvement in precision of the measurement presented here is 
clearly visible. The results confirm that the contribution of the gluon to photon structure is 
significant. The gluon density tends to rise with decreasing x~. The result is compared with 
various parton distributions: LAC1, GRV, SAS1D. They are generally in agreement with the 
data. 



7 Conclusions 

The differential cross-sections da/dpif,, for \rj\ < 1 in the HERA laboratory frame, and da/dr] 
for pt > 2 GeV/c and pt > 3 GeV/c have been measured in photoproduction events with the 
HI detector. The pt spectrum exhibits a high pt tail, larger than in pp collisions at similar 
centre of mass energy. The rj spectra show sensitivity to the parton densities in the photon. 
When charged particles down to pt = 2 GeV/c are included, there is also a significant effect 
from the activity in the underlying event. Charged tracks have been used to extract information 
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Figure 7: The LO gluon distribution in the photon from charged tracks (full circles; < p\ >= 38 
GeV 2 /c 2 ) and jets (open circles; < p\ >= 75 GeV 2 /c 2 ). The inner error bars are statistical, the 
full bars are the statistical and systematic errors added in quadrature. An overall uncertainty 
of 5% from the luminosity measurement and the electron tagger acceptance is not included in 
the errors. The curves are GRV @] (full line), SAS1D |22| (dotted line) and LAC1 f2| (dashed 
line) . 

on the hadronic structure of the photon, by measuring the x 7 distribution. The x 7 distribution 
has been unfolded from the data using for the first time high pt charged tracks to extract the 
LO gluon density of the photon. The gluon density is found to increase with decreasing in 
agreement with an earlier HI analysis using jets to tag and reconstruct the hard sub-process. 
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